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1 Introduction 

Starting from the early seventies, Painleve equations have been playing an increasingly im- 
portant role in mathematical physics, especially in the applications to classical and quantum 
integrable systems and random matrix theory [12], [19], [27], [35], [37]-[40], [42]. A con- 
siderable progress has been made since then in the study of various analytic, asymptotic 
and geometric properties of Painleve transcendents. The interested reader is referred to 
[7, 9, 13] for details and further references. 
The sixth Painleve equation (PVI) 

m » = 1(1 + —!— + —!—) {w f .(l + J- + J-) y + (i.d 

2 \w w-1 w-tJ K ' \t t-1 w-tj v ' 

i 2w(w-i)(w-t) (, in a elt elit-i) {e 2 t -\)t(t-i) \ 

t 2 (t-l) 2 [V 00 2) w 2 (w-1) 2 (w-t) 2 J' 

is on the top of the classification of 2nd order ODEs without movable critical points. 
The latter property means that w(t) is a meromorphic function on the universal cover of 
P 1 \{0, l,oo}. Four complex parameters 6 = (9q, 6t, 6±, 9oo) and two integration constants 
form a six-dimensional PVI parameter space Ai. 
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The most natural mathematical framework for Painleve equations is the theory of 
monodromy preserving deformations. For example, equation (1.1) is associated to rank 2 
system 

dz 

0,t,l) 



1.2 



v z z — t z — 1 J 
with four regular singular points 0, t, l,oo on P . Traceless 2x2 matrices ^4.^ (z/ 

are independent of 2 and have eigenvalues +6 V that coincide with PVI parameters. Also, 
def 

Ao + At + Ai = -Aoc = diag{-6 00 ,6 00 }. 

The fundamental matrix solution <&(z) is multivalued on P 1 \{0, 1, t, oo}, as its ana- 
lytic continuation along non-contractible closed loops produces nontrivial monodromy. Full 
monodromy group is generated by three matrices A^o,t,l £ G = SL{2, C) which correspond 
to the loops 7o,t,i in Fig. 1 (note that MooMiMtMo = 1). Right multiplication of <!> by 
a constant matrix gives another solution, and therefore {.M^} are fixed by (1.2) only up to 
overall conjugation. 

As is well-known, Schlesinger equations of isomonodromic deformation of (1.2) 

dA = [AuM d Ai = [At, Ax] 
dt t ' dt t - 1 

are equivalent to (1.1). Parameter space Ai may be identified with G 3 /G and many ques- 
tions on Painleve VI can be recast in terms of monodromy. This approach, characteristic 
for classical integrable systems in general, turns out to be quite successful. In particular, it 
represents the key element of the solution of PVI connection problem [18], as well as of the 
construction [5, 11] and classification [25] of algebraic solutions. 




Logarithmic derivative of Painleve VI tau function 

o(t) = t(t - 1) — lnr = (t - 1) trAoAt + t tr A t Ai 



(1.3) 



can be expressed in terms of t,w(t) and w'{t) (see e.g. Eq. (2.9) in [26]). It solves a 
nonlinear 2nd order ODE called cr-form of Painleve VI (crPVI), which may be written as 



t{t - l)a' 



-2 det 



261 



to' 



a 



a' + 91 + 6? + 



9 2 

oo 



ta' - a 26f (t - l)a' - a 

\a> + 61 + 01 + 61-91 (t-l)a'-a 26f 



(1.4) 



It is much more natural to work with r(i) than with w(t) for a number of reasons: 
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• First, it is the tau function which typically shows up in the applications of Painleve 
equations, representing gap probabilities in random matrix theory, correlation func- 
tions of Ising model and sine-Gordon field theory at the free-fermion point etc. 

• The notion of tau function extends to the general isomonodromic setting [20] and 
its mathematical meaning is rather clear: it gives the determinant of a Cauchy- 
Riemann operator whose domain consists of multivalued functions with appropriate 
monodromy [34]. 

• Thirdly and (arguably) most importantly, the tau function has an intimate connection 
with quantum field theory. 

The last point was discovered by Sato, Miwa and Jimbo in the first two papers of the 
series [35]. There it was shown that the Riemann-Hilbert problem for rank r linear systems 
with an arbitrary number of regular singularities on P 1 admits a formal solution in terms of 
correlation functions in the theory describing r free massless chiral fermion copies. Besides 
fermions, the correlators involve local fields of another type (below we use the term "mon- 
odromy fields" instead of SMJ's "holonomic") which can be seen in the operator formalism 
as Bogoliubov transformations of the fermion algebra ensuring the required monodromy 
properties. General isomonodromic tau function was originally defined as the correlator of 
monodromy fields. In particular, for Painleve VI it is given by a four-point correlator 

T {t) = (0 Mo (0)0 Mt (t)0 Ml (l)0 Mao (oo)}. (1.5) 

We would like to put isomonodromic deformations into the context of subsequent de- 
velopments in conformal field theory [4, 10, 44]. To our knowledge, no such attempt has 
been made so far except for a short general discussion in [31]. The present paper mainly 
deals with Painleve VI which represents the simplest nontrivial example of isomonodromy 
equations. 

It will be argued below that the relevant chiral CFT has central charge c = 1 and 
monodromy fields in (1.5) are Virasoro primaries with conformal dimensions A v = | tr A„ = 
dl, where v = 0, 1, t, oo. Therefore the structure of the expansion of r(t) near, say, t = 0, is 
strongly constrained by conformal invariance. In fact, the chiral correlator of four primary 
fields for any c is given by 

(to(O)&(*)0i(l)M°°)) = J2 C 0t C ioot Ap ~ A °~ At H^ Ap,q*), (1.6) 

p 

where the sum runs over conformal families appearing in the OPE of 0o and <j>t, A p denotes 
the dimension of the corresponding intermediate primary field (fi p and A = (Ao, At, Ai, Aqo) 

oo 

is the set of external dimensions. Conformal block J-(A, A p ,c;t) = 3~t.(A, A p , c) t 

associated to the channel p is a power series normalized as J~o(A, A p , c) = 1. It is completely 
fixed by conformal symmetry [4]. The structure constants C^ t , C poo combine conformal 
blocks into correlation functions of specific theories and should be obtained from another 
source. CFT correlators usually contain contributions of holomorphic and antiholomorphic 
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conformal blocks subject to the requirement of invariance under the action of braid group on 
the positions of fields. The chiral correlators (1.5), (1.6) are not invariant under this action 
but transform in a natural way, induced by the Hurwitz action on monodromy matrices. 

Direct computation of the coefficients -7-fc(A, A p , c) becomes quite laborious with the 
growth of k. However, very recently this problem was completely solved in the frame- 
work of AGT conjecture [3] relating Liouville CFT and N = 2 4D supersymmetric gauge 
theories. The latter correspondence produces conjectural combinatorial evaluations of con- 
formal blocks, subsequently proven by Alba, Fateev, Litvinov and Tarnopolsky [2]. For 
c = 1, which is the only case of interest for PVI, another derivation was given by Mironov, 
Morozov and Shakirov in [30]. In particular, AGT representation expresses the contribution 
of fixed level of descendants of (fi p to 4-point conformal block ^(A, A p , c; t) in terms of sums 
of simple explicit functions of A, A p and c over bipartitions with a fixed number of boxes. 

Hence, to obtain full expansion of PVI tau function it suffices to determine the dimen- 
sion spectrum of primaries present in the OPEs of monodromy fields and the associated 
structure constants. To formulate the final result, we introduce monodromy exponents 
a = (a ot ,a u ,a i) by 

p^y = tr M^My = 2 cos 271-0-^, /i, v = 0,t, 1. 

Together with 0, these parameters define seven invariant functions on the space A4 of 
monodromy data subject to a relation [18] 

Pot + Pit + Poi + PotPitPoi + Po + Pt + Pi + pL + PoPtPiPoo = (1-7) 

= (POPt + PlPoo)P0t + (PlPt + P0Poo)Plt + (POPI + PtPoo)P01 + 4, 

where p v = tr A4 U = 2cos2-/r0^ [y = 0,t, l,oo). This of course agrees with the dimension 
of Ai, and allows to interpret the triple it as a pair of PVI integration constants. Below we 
assume that 0, a are generic complex numbers verifying Jimbo-Fricke relation (1.7). 

Let Y be the set of all partitions identified with Young diagrams. Given A G Y, we 
write Aj and A^- for the number of boxes in the ith row and the jth column of A, and denote 
by |A| the total number of boxes in A. The quantity h\(i,j) = \'j — i + Aj — j + 1 is called 
the hook length of the box £ A. 

Our main statement is the following 

Claim. Complete expansion of Painleve VI tau function near t = can be written as 

r(t) = const • ^ Cn (0, tr) t^ 0t+n ^- e ^B{6, a ot + n;t). (1.8) 

The function 13(6, a; t) is a power series in t which coincides with the general c = 1 confor- 
mal block J- (0q,9i, 9f , O^Q, a 2 , 1; t) and is explicitly given by 

B(6,a;t) = (l-t) 2e ^ £ (0, <x) i^+W (1.9) 

A l( ueY 
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% + a + i-jf -el) ((e 1 + a + i-j" 2 *2 



~ f — T2 " X 

(ij)eA /i? (i , j) ( A', - i + /ii - j + 1 + 2a 



x 



n 



? t -^ + *-j) 2 -# 2 ) ((^i-^ + ^-i) 2 -^ 



T/ie structure constants {C n (0, (r)} n eZ can be written in terms of Barnes G-function, 
Uee'=± G { 1 + ®t + e6 + ^(ciot + n)) G(l + 0i + + e'(a 04 + n)) 



C n (0,<r) = s 

with s given by 



G(l + 2(a ot + n))G(l-2(a ot + n)) 

(1.11) 



s ±x (cos2n(9 t =F &ot) ~ cos2tt9 ) (cos27r(#i =f cr ot ) - cos 2^^) = (1-12) 
= (cos 2n9 t cos 2ir6i + cos 27r#o cos 2^6*00 ± i sin 27rcrot cos 2iraoi) — 
— (cos 2-kOq cos 27r^i + cos 2ir9t cos 2^9^ =p z sin 2iraot cos 2nait) e =l:27rio "o t . 



Analogous expansions of r(t) at t = l,oo are obtained from (1.8)-(1.12) by applying a 
suitable transformation of parameters. 



Painleve transcendents are generally believed to be rather complicated special functions. 
In our opinion, this reputation is somewhat undeserved. Painleve VI, for instance, enjoys 
most of the basic properties of the Gauss hypergeometric equation. In particular, it has 
many elementary solutions (see [25] and references therein), Backlund transformations [33], 
quadratic Landen-type transformations [21], and its connection problem is solved [18]. The 
present work adds to this list a link to representation theory of the Virasoro algebra and 
the series representation of PVI solutions, which can also be seen as an efficient tool for 
their numerical computation. 

The outline of the paper is as follows. The next section starts with a brief survey of the 
isomonodromy problem. After exhibiting global conformal symmetry of the tau function, 
in Subsection 2.3 we introduce monodromy fields and explain how various mathematical 
objects of the theory of monodromy preserving deformations can be written in terms of 
correlation functions in 2D CFT. Subsection 2.4 deals more specifically with Painleve VI. 
Here we present the arguments leading to (1.8)— (1.12), and discuss analytic continuation 
of PVI solutions and their Backlund transformations from CFT point of view. Section 3 is 
devoted to direct analytic verification of the above expansion. This task is further pursued 
in Section 4, where the known special PVI solutions are examined from the field-theoretic 
perspective. We conclude with a list of open questions and directions for future work. 
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2 CFT approach to isomonodromy 



2.1 General 

It is instructive to start with the general case of rank N linear system with n regular singular 
points a = {ai, . . . , a n } on P 1 . Instead of (1.2) one has 

n A 

The absence of singularity at infinity implies that N x N constant matrices {A u } satisfy the 
constraint Ylu=i^ = They are assumed to be diagonalizable so that A v = Q v T v Qy 
with some T v = diag {A^i, . . . , A^jv}. The fundamental solution will be normalized by 
$(^o) = IjV- It is useful to introduce the matrix 

J(z) = S" 1 ^ = $~ 1 ^)$. 

The coefficients of the Taylor series of &(z) around z = zq can be expressed in terms of J 
and its derivatives. In particular, 



( — ) 2 

${z -^z ) = l N + J (z ) (z - zq) + {J 2 {z ) + dj(zo)) 1 . ° ; + . . . (2.2) 



Near the singular points, the fundamental solution has the following expansions (under 
additional non-resonancy assumption \ v a — \ v k £ % for j ^ k): 

§{z^a v ) = Q v {z)(z-a v ) Tv C v . (2.3) 

Here Q v {z) is holomorphic and invertible in a neighborhood of z = a u and satisfies Q u {a v ) = 
Q v . The connection matrix C v is independent of z. Counterclockwise continuation of $(z) 
around a v leads to monodromy matrix Ai u = C~ 1 e 2nz ^C 1/ . 

Let us now vary the positions of singularities and normalization point, simultaneously 
evolving A^s in such a way that the monodromy is preserved. A classical result translates 
this requirement into a system of PDEs 

8^9=-?^-^-*, (2.4) 
zq - a u z - a u 

d Z0 <S>= -A(z )$. (2.5) 

It is important to note that the matrix J{z) remains invariant under isomonodromic vari- 
ation of zq. Schlesinger deformation equations are obtained as compatibility conditions of 
(2.1), (2.4) and (2.5). Explicitly, 

d a A v = - , n ^ u, (2.6) 

z - a M a M - a v 

^x = -V^^, daa A v = -Y)^A. (2.7) 
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Lax form of the Schlesinger system (2.6)-(2.7) implies that the eigenvalues of A^s are 
conserved under deformation. This is of course expected due to the obvious relation between 
the spectra of A„s and monodromy matrices. 

Isomonodromic tau function r(a) is defined by 

rflnr = tr A^A U din (a^ — a u ) . (2.8) 

It is a nontrivial consequence of the deformation equations that the 1-form in the r.h.s. is 
closed. To show that it does not depend on zo, one can rewrite (2.8) as 

r = £ tr^A = 1 tlJ 2 {z y (2 9) 

Here the first equality follows from (2.8) and the second one from the fact that 3{z) is 
conjugate to A(z). 

Finally, let us decompose all ^'s into the sum of scalar and traceless part as A v = 

tr A - 

ljv + A p . It can then be easily checked that 



v v 

ta(o) = J] (a„ - aS^^ r A {a). (2.12) 

jl<V 

This allows to assume without any loss of generality that A^s are traceless, but we delib- 
erately postpone the imposition of this condition. 

2.2 Global conformal symmetry 

az + p 

Fractional linear maps f(z) = form the automorphism group of the Riemann sphere. 

72; + 

It is clear that under the action of these transformations on z, zq and a the quantities $(z) 
and {A u } transform as functions and J{z) as a vector field. Our task in this subsection is 
to understand the effect of global conformal mappings on the tau function. 

Let us first compute r(a) explicitly in the case n = 3. Since Ai + A2 + -4-3 = 0, the 
coefficients tr A^Av in (2.8) are conserved quantities. Hence, integrating (2.8), we find 

r( ai ,a 2 , 03) = const • - a 2 ) A3 " A2 - Al (ai - a3) A2 - Al - A3 (a 2 - aa)^ 2 "^ , 

where A u = ^trAl with v = 1,2,3. One recognizes here the general expression for the 
three-point correlation function of (quasi) primary fields of dimensions Ai t 2,3 in the two- 
dimensional conformal field theory. 

Given the above example, it is natural to assume that for general n the tau function 
transforms as the n-point function of primaries with appropriate dimensions: 



T(f(a)) = H[f>(a u )]- A » r(a). (2.13) 
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To prove the last formula, it is sufficient to consider infinitesimal transformations generated 
by the vector field Z = (A + Bz + Cz 2 ) d Zl which amounts to checking three differential 
constraints: 

^d a „ lnr = 0, 
5^(a„d a „ lnr + A„) = 0, 
^(a*A, lnT + 2A u a u ) =0. 

V 

These relations can indeed be straightforwardly demonstrated using the first equality in 
(2.9) and the condition ^j V A v = 0. 

2.3 Field content 

The fundamental solution $ is completely fixed by its monodromy, normalization and sin- 
gular behaviour (2.3). The last point is particularly important: indeed, adding an integer 
to any diagonal element of T u modifies the asymptotics of as z — > a u without changing 
monodromy matrices. Therefore, in what follows, the notion of monodromy will include 
not only the set of Mu s but also the choice of their logarithm branches C u = C~ l T u C u . 

Let us now try to construct a formal QFT solution of the isomonodromic deformation 
problem, extending the ideas of [31], [35]/_77. The starting point will be the following 
ansatz for <£: 

w - (, - *r {0ci { T ' v °f * (2>) . (2.i4) 

{Oc 1 {a 1 )...Oc n {a n )) 

Here it is assumed that {OcJ\, {<fik} ar e primary fields in a 2D CFT characterized by 

some central charge c. Further, we want the OPEs of <^'s with </?'s to contain the identity 
operator. This forces them to have equal dimensions, to be denoted by A. Normalization 
of these fields is fixed by the normalization of <£; the leading OPE term should be equal to 

^Pj{zo)^k{z) ~ (z - z )~ 2A S jk . (2.15) 

Since $(z) may be represented by an entire series near z = zq, the dimensions of all other 
primaries appearing in this OPE should be given by strictly positive integers. Monodromy 
fields {Oc v } are defined by the condition that their complete OPEs with {ipk} have the 
form 



Cv (au)^Pk(z) =J2({z- a v ) c ") {z 



a 



jk ■ 

j=i e=o 



V) J 



where {Oc„,j,i} are some local fields. In particular, the raw vector (ip\ . . . <p n ) should be 
multiplied by M v when continued around Oc v ■ If one succeeds in finding a set of fields with 
all mentioned properties, the correlator ratio (2.14) will automatically give the solution of 
the linear system (2.1). 

The definition (2.8) of the tau function arises very naturally in the CFT framework. 
To illustrate this, let us compute two more orders in the OPE (2.15). The identity field 
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has no level 1 descendants, therefore the leading correction is given by a new primary field 
Jjk of dimension 1. The next-to- leading order correction comes from three sources: i) 
nonvanishing level 2 descendant of the identity operator given by the energy-momentum 
tensor T, ii) level 1 current descendant dJj k and iii) new primaries of dimension 2 which 
can be combined into a single field Sjk- Thus 

Vj{zo)<Pk{z) = (z - z y 2A 5 jk + Jjk(zo) (z - Zq) + (2.16) 

+ (^ T(z )S jk + (dJ jh )(z ) + S jk (z )) {Z ~ 2 Z ° )2 +0((z- z f) 

We will make a further assumption of tracelessness of S, which is essentially motivated 
by the examples considered below. Now, substituting (2.16) into (2.14) and matching the 
result with (2.2), one finds that 

(0 Cl (ai)---Oc n MJ(z)) 

tr^^^ ^f^f^ 4 ^. (2.18) 
Standard CFT arguments allow to rewrite the r.h.s. of the last formula as 

(0 Cl (a 1 )...Oc n (a n )T(z)) _ A \ K 1 . f , , S\ 

— Tr\ — 1 — ^ rs — 1 — vi — 2^ \ 7 \2 H da - in ( a i) • • • u C n K a n)) > , 

(0 Cl {ai) ■ ■ ■°c n { a n)) ~i\^ z ~ ay ' z ~ av J 

where Aj, denotes conformal dimension of Oc v - Comparison of (2.18) with (2.9) then shows 
that the tau function can be identified with a power of the correlator of monodromy fields, 

2NA 

r(a) = (O jCl (a 1 )...Oc n (a n ))— . (2.19) 

In what follows, we will be exclusively interested in the case when 

c = 2iVA. (2.20) 

Such a condition implies, in particular, that the dimensions A u of monodromy fields coincide 
with the quantities A„ = | tr Al from the previous subsection. 

One possible realization of the above conditions is provided by the theory of iV free 
complex fermions. Its central charge c = N agrees with the conformal dimension A = ^ of 
fermionic fields {ijjj}, {ipk} which play the role of (p's and (p's. The currents are by definition 
given by Jjk = (ipj ipk), while the energy-momentum tensor T and the fields {Sjk} may be 
expressed as 

k 

2 

Sjk = \i>3 dipk) - {dipj tp k ) - —T5j k . 
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To represent monodromy fields, recall the usual bosononization formulas 




2 

k 

where {(j>k}k=l,...,N ar e free complex bosonic fields with the propagator {cj>k(w)(j)k(z)) ~ 
— ln(z — w). Also note that for C £ GL(N,C), monodromy matrices for the linearly 
transformed fermions 

k k 

are obtained from Ai^s by conjugation by C. In particular, setting C = C u , one obtains 
fermions {^i^}, {V^i with diagonal monodromy around a v . Denote by {</>^' > } bosonic 
fields associated to this "diagonal" fermionic basis, then monodromy field Oc v can be written 
as 

We thus need to deal with n different bosonization schemes of the same theory, each of them 
being adapted for representing one of the monodromy fields. The corresponding iV-tuples 
of bosons are related by complicated nonlocal transformations. 

The formulas (2.10)-(2.12) are a signature of the well-known decomposition of fermionic 
CFT into the direct sum u(l) © su(N)i of two WZW theories. Fermion and monodromy 
fields are given by products of fields from the two summands: 

Cv = :e ^ 00 :®O l . 

Bosonic field <pQ in the u(l) factors is expressed in terms of fields introduced before as 

<f)Q = — = "^2^=1 §k- The fields {<fik}i {<fik} an d {O £ } live in the su(N)i WZW theory and 

can be formally written as ordered exponentials of integrated linear combinations of su(N)i- 

currents. It should be emphasized that they are Virasoro primaries but not necessarily 

N-l 

WZW primaries. The fields {tfk} and {(f>k\ have the same dimension A = , in 

accordance with the central charge c^nyv = iV — 1 [23]. The dimension of O £ is equal to 

tr A. 

2 tvA„, where as above, A v = A u ljv stands for the traceless part of A v . 

Now it becomes clear that imposing the tracelessness of A(z) corresponds to factoring 
out the u(l) piece from the fermionic theory. This innocently looking procedure is in fact 
crucial, as it drastically reduces the number of primary fields in the OPEs and thus makes 
the computation of correlation functions much more efficient as compared to fermionic 
realization. Therefore, in what follows we set tvA(z) = 0, remove the hats from A^s, C u 's, 
0's and ip's to lighten the notation, and interpret the isomonodromic tau function as a 
correlation function of primaries with dimensions A u in a CFT with c = N — 1. 
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We close this subsection with an example of application of field-theoretic machinery 
in the case N = 2. It is somewhat distinguished from CFT point of view, since for c = 1 
the dimension A = | of </?'s and tp's corresponds to level 2 degenerate states, and the 
dimension 1 of {Jjk} is degenerate at level 3. Hence the correlation functions 

Vjk = • • • Cc n (a n )ifj(zo)ip k (z)}, 

Qjk = {Oc l (a 1 )...Oc n (a n )J jk (z)), 

have to satisfy linear PDEs of order 2 and 3, fixed by Virasoro symmetry. This results into 
the following statement (cf observations made in [32]): 

Proposition 1. Under assumption tr A(z) = 0, the matrices 

V = (z-2 )-5T$, Q = T <$>- 1 d z <&, 

satisfy the differential equations 

d zz V=\^—d Z0 + 1 + W— L_ d av + A \ 2 )\v, 
\z-zo 4:(z-z ) 2 ^\z-a v ( z - au ) 2 Ji 



[ v\ z - a v {z-a u ) J e -f\\z-a v ) (z-a v )J] 



Proof. Straightforward but tedious calculation using the relations (2.1), (2.4), (2.5), (2.9) 
and the identity A 2 = \ tr A 2 I2 verified by any traceless 2x2 matrix A. □ 



2.4 Painleve VI 

Recall that global conformal symmetry allows to fix the positions of three singular points. 
Painleve VI equation corresponds to setting N = 2, n = 4 and sending these three points to 
0, 1 and 00. The remaining singular point, z = t, represents the cross-ratio of singularities, 
which is preserved by Mobius transformations. 

For v = 0, t, l,oo, let us denote by ±9 U the eigenvalues of A u . Preceding arguments 
show that PVI tau function r{t) defined by (1.3) is nothing but the four-point correlator 
of monodromy fields, 

r(t) = (O C0 (0)O Ct (t)O Cl (l)O Coo (^)), (2.21) 

and that these fields are Virasoro primaries with dimensions A u = 9 2 in a c = 1 confor- 
mal field theory. The field at infinity should be understood according to the usual CFT 
prescription 

(...0{oo)) d = f lim R 2A °(...0(R)). 

It is clear that auxiliary fields {ifk} should have monodromy AifMo around all fields 
in the OPE of Oc and 0£ t . Let e ±2m(T0t denote the eigenvalues of MtM-o and Cot be its 
diagonalizing transformation. Since o~Qt is defined only up to an integer, it is natural to 
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expect that the set of primaries present in the OPE of Oc and Oc t consists of an infinite 

number of monodromy fields O ( n ) with n G Z and 

'-ot 

C tt = C 0t ( ""n " n I C °*' 

i.e. of all possible monodromy fields associated to the monodromy matrix .Mi-Mo- Taking 
into account that conformal dimension of O (») is equal to (o"ot + n) 2 , the first part of 
our main statement (formulas (1.8)-(1.10)) now follows from the general formula (1.6) and 
AGT combinatorial representations of conformal blocks [3]. 

The structure constants C n (6,cr) of the expansion (1.8) can be determined from the 
so-called Jimbo asymptotic formula [18], expressing the asymptotics of PVI tau function 
as t —7- in terms of monodromy. In fact we have already obtained the "easier half" of this 
formula. E.g. if — g < Reuot < i, then (1.8) implies that the leading behaviour of r(t) is 
given by 

r(t -> 0) ~ const • t^ 9 '^ . 

Subleading asymptotics, fixing the second PVI integration constant, can be rewritten in 
the form of a recursion relation on the coefficients C n (0,cr). Namely, 

C n ±i _T 2 (1 T 2(a 0t + „)) ^ r (1 + e9 + 6 t ± (a t + n)) V (1 + eO^ + 01 ± (ff M + n)) 

■ X 



n 



C n T 2 (1 ± 2(a ot + n)) ^ r (1 + e6 + 0tT ((Tot + n)) r (1 + + 0! t (^ + n)) 

x (gg ~ (gt =F (got + n)) 2 ) (C ~ C?i =F (fot + r^)) 2 ) ±1 
4(a 0i + n) 2 (l±2( ( To i + n)) 2 Sj ' 

where s is defined by (1.12). This relation can be easily solved in terms of Barnes functions, 
with the answer given by (1.11). It is interesting to note that, up to a common multiplier 
and appropriately symmetrized s n factors, C n 's essentially coincide with the chiral parts 
[36] of the corresponding structure constants in the time-like Liouville theory [16, 43]. 

Remark 2. The structure constants (1.11) can not be completely factorized into the prod- 
ucts of three-point functions due to the presence of the parameter s. This is an artifact of 
non-trivial braid group action on the correlation functions of monodromy fields. 

To illustrate what we have in mind, consider the analytic continuation of r(t) along 
a closed counterclockwise contour around the branch point t = 0. In general, such a 
continuation induces an action of the 3-braid group (more precisely, of the modular group 
r(2)) on monodromy [11]. In the case at hand, new monodromy matrices are given by 

M = M t M Mi\ M' t = (MtMo)Mt(M t M )- 1 , M[ = Mi, 
so that a' ot = <7ot and 

Poi =PoPi + PtPoc - Poi - PotPlt, (2-22) 
Pit = PWt + PoPoo - Pit - Potp'oi ■ (2.23) 
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Therefore, the change of the branch of r(t) is encoded in the change of the structure con- 
stants. On the other hand, one can perform analytic continuation directly in the expansion 
(1.8). Up to an irrelevant overall factor, this amounts to multiplication of C n (0,cr) by 
e 47rm<T(tt . Since both results should coincide, the structure constants have to satisfy the 
functional relation 

c n (e,<j') = K-e 4 ™ n ^c n {e,cT), 

where k is independent of n. The factor s n in (1.11) is a minimal solution of this relation, 
as for a' defined by (2.22)-(2.23) one has s(0,tr') = e 4 ™ ot s(0, a). 

Remark 3. Let us denote by dim A the number of standard Young tableaux of shape A G Y. 
It coincides with the dimension of the irreducible representation of symmetric group S\\\ 
associated to A. Also write d\ for the number of diagonal boxes in A and introduce the 
Frobenius coordinates 

Pi =\~i, Qi = A- - i, i = l,..., d x , 

which give the number of boxes to the right and above the ith diagonal box. It is well 
known/easy to show that 

dim A 1 1 

det 



Ui^exhxihj) n^ir(^ A + i)r(^ + i) 



Pi + Qj + 1 



Y[ (i-j+ z )(i-j+3>) = (zzr) d *i[T 

(m)6A i=l 



i,j=l,...,dx 

p^ + l + z, qi+1- z,p^ + 1 + z', q£ + 1 - z' 
1 + z, 1 - z, 1 + z', 1 - z' 



We now recognize in (1.10) typical pieces of z-measures on partitions [6]. It would be 
nice to understand this coincidence conceptually with the purpose to sum up the series for 
B(0,<r; t) and r(t). 

Remark 4. Painleve VI equation has a large group of hidden symmetries of affine Weyl 
type [33]. Almost all of them are manifest in the conformal expansion (1.8). For instance, 
the change of sign of any of parameters 6 has no effect on the tau function, since conformal 
blocks depend only on the dimensions A u = 9% and the ratios of structure constants (1.11) 
also remain invariant. 

Conformal block symmetry (0o,0t) <->■ {0oo,0i) and its counterparts for expansions at 
t = 1, oo yield further simple transformations of r(i). Another, less trivial symmetry that 

can be found by inspection of (1.10) shifts the values of all 6 by 5 = - — 1 1 — . 
Additional transformations come from crossing symmetry. In contrast to the previous ones, 
they also act on t by fractional linear transformations exchanging 0, 1 and oo. 

The action of generators of the above transformations on the parameters 6, conformal 
blocks and tau function expansions is recorded in the following table: 
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#0 


0t 


01 


0oo 


B(G,a;t) 


r(t) 


*0 


— On 
u 


0* 

u t 


u l 




u oo 


B(Q rr-t) 


rff) 


St 


Oo 


-0t 


01 


0oo 


B(0,a;t) 


r(t) 


Si 


Oo 


0t 


-0i 


0oo 


B(0,a;t) 


r(t) 


Soo 


00 


0t 


0i 


— 0oo 


B{0,a;t) 


/ i \ 

r(t) 




On - 5 


0+ -5 


Oi - 5 


"OO u 


(I _ t) 5ltS B(0 a- 1) 


t 5 0t S (\ — t) SltS T(t) 


rot 


a 

Poo 


a 
"l 


a 
vt 


Oo 


d(U, a , t) 


I T[tj 


' It 


0+ 


On 
a 


a 

"oo 


0, 

"1 




{1 I) T{Z) 


roi 


01 


0oo 


#0 


0t 


(l-t) A "£(0,a;t) 




9oi 


01 


0t 


00 


0oo 




r(l-t) 


90oo 


0<x> 


0t 


01 


00 






<7loo 


00 


0t 


0<x> 


0i 


(1 — i) A o-A t -A CTj g {q^ a . 


( i-r 2A ^(^) 



Here we have introduced the notation 



Sot — 0o + 0t — 0i — 0oo, — Ao + A t — Ai — Aqo, 

Su = Oi + 0t-0o- 0oo, A w = Ai + A t - A - Aoo, 

and A CT = a 2 . Transformed tau functions in the last column are (in some cases) defined up 
to constant factors which depend on the choice of normalization of the structure constants. 

To complete the picture, it remains to understand the QFT meaning of an elementary 
Schlesinger transformation, e.g. the one shifting Oq and Oi by i. It may be expected that 
this symmetry arises from the fusion of auxiliary fields <fj and (fk with monodromy fields 
Oc and Oc± in the correlator representation (2.14) of the fundamental solution <3?. Indeed, 
<^'s and (£>'s are degenerate at level 2, and therefore their OPEs with monodromy field of 
dimension 2 can only contain two conformal families generated by monodromy fields with 
dimensions (0 db ^) 2 . 



3 Painleve VI recurrence 

Painleve VI tau function expansion (1.8)— (1.11) is, in the strict mathematical sense, a 

conjecture. On the other hand, r(t) is completely fixed by its leading asymptotics as t — > 0. 

Ci(0 <x) 

Once the quantity not in the power-law exponent and the amplitude ratio 1 — - are 

G {0,(T) 

found from Jimbo's formula, the rest of the series can, at least in principle, be recursively 
reconstructed order by order from uPVI equation and checked against CFT predictions. 
Below we describe technical details of this procedure and derive several first terms of the 
conformal expansion. 

The equation (1.4) gives a quadrilinear 3rd order ODE for the tau function itself. 
Observe, however, that if we differentiate (1.4) one more time with respect to t, the resulting 
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equation will be divisible by cr"(t) and in the end turns out to be bilinear in r(t). It is also 
convenient to work with a slightly modified function rf(t) = £ Ao+A *(l — i) At+Al r(i), which 
partially takes into account the behaviour of r(t) as t — > 0, 1. It satisfies the equation 

t(t-l)r] 2 \ -8[t(t-l)(V) 2 J -4(2A-l) W " + 8(A-2)(7/) 2 }+ (3.1) 
+ m 3 {t - l) 3 (t]") 2 + 4t(t - 1)7777" - 4(,4 - - l)(2f - 1)7777' + A{Bt + C)?? 2 = 0, 



where A, B, C are given by 

A= Ao + At + Ai + Aoo, 
P = (A - Ai) (Aoo - A*) , 
C = (A - A t ) (Ai - Aoo) , 

and the fourth PVI parameter is killed by differentiation. 
Now substitute into (3.1) the ansatz 



00 

2 ^ (n). k 



(3.2) 



fc=0 



(n) 

normalized by ?7q = 1 for all n £ Z. Picking up the coefficients of different powers of t in 
the result, one obtains an overdetermined system of equations for {C n }, {77^}. Namely, 
for any I E Z and any L £ Z>o we have 



m,nd1i 
m+n=£ 



V r r V rv fcfc 'r7 (m) r7 (n) 4- V R kk ' ri {m) ri^ + 

/ j ^m^n / J ^mn'lk 'Ik' ' / i ^mn'lk Ik' ' 



kk'{m) (n) 



k,k'>0 
k+k' +m 2 +n 2 =L+l 



k,k'>0 
k+k'+m 2 +n 2 =L 



+ E E W»£ 



kk'(m) (n) 



k,k'>0 
k+k' +m 2 +n 2 =L-l 



k,k'>0 
k+k'+m?+n 2 =L-2 



0. 



(3.3) 



with somewhat cumbersome but completely explicit polynomial coefficients: 



kk' 



kk' 



linn 



\kk' 



kk' 



kk 2P - f kk 4- 1 



= 3 e 



.fcfe' 



2 (P - 1 + 2,4) e£i - P (P - 2A) + 4C, 



3 e 



.fcfc' 



2(P-4A)e*£ + (P-1)(P-1-2A) + 4B, 



kk' 



2P + e 



A- A:' 



3-4A 



k-k' + {rn-n){l + 2a)\ 



P = 2a 2 + 2(j^ + L. 



Let us write down the recurrence relations (3.3) and their consequences for some L 
and I. 
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L = 0, £ = : from (3.3) it follows that 



which in turn gives 



rtf> = -(A t + Ai) + (A - " A ° + A ' )(Act " A - + Ai: 



2A a 

with A CT = o" 2 as above. This reproduces n = 0, level 1 descendant contribution to the 
expansion (1.8). The corresponding term was already found in [18]. 

L = 1, £ = ±1: in this case, (3.3) implies that 

1,0 . 0,1 \ (±1) . / 0,1 . 1,0 \ (0) . / o0,Q . o0,0 \ n 

"±1,0 + «o!±i )m + ^±1,0 + a o!±i ) m + (p±i,o + Po,±i ) = °> 

and we obtain 

(±i) , A . A v . (A CT±1 -A + A t )(A (T±1 -A oo + A 1 ) 
?/i --(Ai + Ai)H 



2A, 



±i 



The latter expression corresponds to level 1 descendants with n = ±1. It coincides with 
with a replaced by cr±l, but this is not surprising: the same should be correct for any 
ngZ provided the conjectured periodicity of powers in (3.2) holds true. 

L = 1, £ = 0: here we find 



r 2 

u 



a o°o + a ofi) ffl + «o!o (^P) + (A),o + ^o,o) »7i 0) + 7o,o 



+ (3-4) 



+CiC_i fa^° ! + = 0, 

which gives level 2 descendant contribution with n = 0: 

,f = - (A ' + Al)( t + Al + 1) - (A. + A^fh 
(Ag - A + A t )(A g - A + A t + 1)(A CT - Aqq + Ai)(A g - A^ + Ai + 1) 

4A ff (l + 2A a ) 

(1 + 2A(j) (A + A t + A g (A g -l|-3(Ao-AO^ j^ Aoo + Al + ^(A^-lKA.-AO' 

+ 2(1-4A CT ) 2 

To obtain the last formula, one should use, in addition to (3.4) and the coefficients found 
above, the relation 

dCLi (( A ° + A * " A ^) 2 ~ 4A o A t) (( A oo + A x - A CT ) 2 - 4AOOA1 
C~o 16A 2 (1-4A CT ) 2 ' 

This piece of initial data disappears after the above differentiation but can be determined 
from the quadrilinear form of erPVT. 
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It is straightforward to compute more terms in (3.2) using computer algebra. The 
procedure works as follows. For a given L, one should start with maximal \£\ producing a 
nontrivial relation (3.3), and then repeatedly decrease \£\ by 1. When all possibilities are 
exhausted, increase L by 1. The coefficients determined at the first steps of this iterative 
procedure are listed in the table below. Empty entries correspond to the relations satisfied 
automatically or to no relations at all. 








1 


2 


3 


4 


5 


6 


7 


8 


9 





(0) 

Vi 


(0) 
V 2 


(0) 

V3 


(0) 

vl 


(0) 


(0) 

Ve 


(0) 


(0) 

Vs 


(0) 
Vg 


(0) 

Via 


±1 




(±1) 

Vi 


(±i) 


(±i) 
m 


(±i) 

v\ 


>T 


(±1) 


(±i) 

v) 


(±1) 
Vg 


(±i) 
Vg 


±2 










(±2) 

Vi 




(±2) 

m 


(±2) 

v\ 


(±2) 
V5 


(±2) 

m 


±3 


















C±3 


(±3) 

Vi 



In this way, we have successfully checked the expansion (1.8)— (1.11) for n = 0, ±1, ±2, 
±3 going up to level 10 in descendants. To give the reader an idea of the computational 
complexity, we note that there are nearly 500 bipartitions of size 10, each of them producing 
a rational function of 6, aot in the corresponding expansion coefficient. 

4 Conformal blocks and special PVI solutions 

4.1 Riccati solutions 

These PVI solutions appear when the monodromy of (1.2) is equivalent to an upper trian- 
gular one. Parameters 6 can be Backlund transformed to satisfy Oq + 0t + 0\ + Ooo = 0. The 
initial conditions are also constrained and can be chosen as <x = (0q + Of, 9\ + Of, 0q + 0\). 

This results into a one-parameter family of PVI transcendents w(t) that may be written 
in terms of Gauss hypergeometric functions, see e.g. Proposition 49 in [25] for explicit 
formulas. The relevant tau function, however, is extremely simple: 

r(t) = const 

We recognize in the r.h.s. of this relation the four-point correlator (Ve (0)Vg t (t)Ve 1 (1) Vg^ (oo)) 
of the chiral vertex operators V$(z) = : e*^ 6 ^ 2 ) : made of free massless bosons. The contri- 
bution of conformal blocks 13(0, crot + Ti; t) with n ^ to the expansion (1.8) is annihilated 
by the vanishing structure constants. 

4.2 Chazy solutions 

As is well-known [4, 10, 44] and already mentioned in Section 2, the presence of degenerate 
states in the Virasoro module generated by a primary field <ft leads to linear differential 
equations for the correlation functions involving this field. Possible conformal dimensions 
of such (fi's are determined by the zeros of Kac determinant. They are labeled by two integers 
r,s £ Z>x and, for c = 1, are explicitly given by A^ rs = ^ r ~^' . In the simplest nontrivial 
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case r = 2, s = 1 the four-point correlator containing (j>2,i(z) and three arbitrary primaries 
solves a linear 2nd order ODE that can be reduced to the hypergeometric equation. 

Thus if one of the monodromy matrices -Mo,t,i,oo is equal to —1, there exists a PVI tau 
function given by a linear combination of two solutions of the latter equation. For example, 
for 6oo = ±2 (i.e. Aqo = |) the general solution of 



r" + 2 



2+02+1 02_02 + 02 + ! 



+ 



t - 1 



+ 



i(t-l) 



e=± 



1\2 



9 + e x ) 2 + 0f - \ + ^ / flg+g2-( e g 1+ i)2 + 0? + 2 - (, 0, , ! k ) 



t- 1 



T + 



T = 



satisfies uPVI. Conformal expansion of this tau function at, say, t = is determined by 
only two channels with dimensions (#i ± ^) • More precisely, one has 



r(t) = Y,c^Hf-^ B (e o ,e t ,0i, Idh. + ht), 



(4.1) 



where C± are arbitrary constants. The relevant conformal blocks can be written as 



B (<h, 9t,e h 1 -,9 1 + 1 -;t)=(l- t)<«*&-*-* 2F X 



0Q + St + 01 + 2 J ^0 — Qt + #1 + 2 
1 + 20! 



The tau function (4.1) gives another known class of special function solutions of Painleve VI, 
the so-called generalized Chazy solutions (cf Lemma 33 in [29]). 



4.3 Algebraic solutions 

Algebraic PVI solutions correspond to finite orbits of the braid/modular group action on 
monodromy of the associated linear system [11]. All such solutions have been classified 
in [25]. It turns out that there are 45 exceptional equivalence classes with fixed rational 
parameters 9, a and three families continuously depending on some of them. 

In the exceptional cases, the contributions of different conformal blocks overlap in the 
tau function expansion, which makes difficult their identification. On the other hand, this 
task is rather straightforward for continuous families, see examples given below. An infinite 
number of other explicit examples of c = 1 conformal blocks can be generated by Backlund 
transformations. 

Example 5. For 6 = (a, a, b, \ — 6), a = (2a, |, j) there is a solution 



r(t) = const • t 2a (1-t) 



-a 2 -b 2 - 



B(a,o,6,i-6,2a;t) 

arising from the contribution of a single conformal block (Solution II in [25]). 



-18- 



Example 6. For = (2a, a, a, g), tr = (3a, |, |) there is a solution that can be parame- 
trized as follows (Solution III in [25]): 



2_ 1 



r(t(s)) = const • (a - 2) 4a (a - iy a ~^(a + l)- iUcl +s( s + 2) 
ifal " (g + ^ " 2) 

Consider the interval s E (2, oo) which maps to t E (0, 1). Again only one block contributes 
to the tau function expansion on the corresponding branch at t = 0. It is explicitly given 
by 

/ i i \ -18a 2 + | / . \ 14a 2 -| 

B(2a,a,a,i,3„; t ) - <„- l)"""* (*±i) (*±?) 



with = - — — — — - ' . Note that i = is not really a branch point of St. 
(l-t)« 



13 



Example 7. Taking = (a, a, a, |), <r = (2a, g, |) one has a solution (corresponding to 
Solution IV in [25]) 



16 . 



r(t(s)) = const ■ a"** + I3(s-l)-~ +H(s+1)^ (2 - »)■*■ (2s - 1) 
w 2a - 1 

Under this parametrization, the interval s E (1,2) is mapped to t E (0,1). Tau function 
asymptotics near t = on the relevant branch yields the conformal block 

1 „„.A-r.. n-6« a +i (st\-^ 2 +h fs t + l\ 2a2 (2s t -l^ 7 " J - J 



B (a,a,a,-,2a;t) =(s t -l 



where 



s t = \ U+n 4 + y3-n2+ 2 — -). „, y'l -(17(1. -/)): 



4.4 Picard solutions 

The remainder of this section deals with Painleve VI solutions of Picard type. Here the 
parameters are chosen to be 0pi car d = (i ; j, j, !)■ Up to Backlund transformations, this is 
the only case where the general two-parameter solution of PVI is available [14, 28]. Quite 
remarkably, for precisely these 6 there exists an explicit formula for c = 1 conformal block, 
found by Zamolodchikov [41]. We now briefly explain the relation between the two subjects. 

First recall that complex torus can be seen as a two-sheeted covering of the four- 
punctured sphere, its period ratio being determined by the cross-ratio of the punctures. 
This suggests an elliptic parametrization of PVI variable t: 

inr iK'(t) 

q=e r= w 

where K(t), K'(t) are complete elliptic integrals of the first kind: 
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In this parametrization, conformal block with external dimensions corresponding to 
0Picard is given by (see Eq. (2.28) in [41]) 



B (#Picard> C; t) 

(l-t)i0 3 (O|r) 

Here and below $3(z|r) = ^2 e mn2T + 2mz denotes Jacobi theta function. On the other hand, 

Barnes function duplication identity reduces the formula (1.11) for the structure constant 
to 



(Cn, 



7T 



Picard 



^ (I) 



(_ s ) n 2~ 4 ( cr °* +n ) 2 



COS 7T(Toj 

The expression for s also simplifies drastically. Jimbo-Fricke relation (1.7) implies that 
for fixed aot, &it the two possible values of cos27rcxoi are given by — cos27r(<7ot =F ffit) and 
therefore s = — e ±2mcFu . 

Summation over conformal families in (1.8) now gives theta function series so that we 
get an explicit expression 



TPicardW = Const 



At 



#3 (omttt ± cri t 7r|r) 



3 (O|t) 



Straightforward check (of the literature [22]) shows that this function indeed satisfies PVI 
with Picard parameters, i.e. the expansion (1.8) is complete. 

The above observations on the correspondence between special solutions of Painleve VI 
equation and c = 1 conformal blocks are collected in the following table: 



Painleve VI solutions 


c = 1 conformal blocks 


Riccati 


vertex operators 


Chazy 


singular vectors 


Picard 


Ashkin- Teller 


algebraic 


? 



5 Further questions 

To summarize, we have obtained a complete series expansion of Painleve VI tau function 
near the singular points t = 0, 1, oo. This series involves summation over conformal families 
arising in the OPEs of monodromy fields and labeled by n E Z, as well as a double sum over 
Young diagrams which encode the contribution of Virasoro descendants. The field theory 
meaning of r(t) now becomes clear: it is a generating function of c = 1 conformal 
blocks. 

This opens a way to intriguing applications of AGT correspondence in the theory of 
monodromy preserving deformations. To increase the number of singular points of the 
linear system (1.2), one should deal with AGT representation for the n-point conformal 
block with n > 4. Increasing rank leads to other integer values of central charge, and it 
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can be expected that intermediate conformal families will be labeled by (N — l)-tuples of 
integers. Another interesting option is the study of isomonodromic tau functions on higher 
genus Riemann surfaces. 

Conversely, isomonodromy problems may provide useful insight into CFT. For instance, 
it would be interesting to understand if conformal blocks associated to exceptional algebraic 
Painleve VI solutions can be computed in explicit form and identify the corresponding 
theories, orbifold CFTs [15] being the most natural candidates. One may also try to generate 
new explicit examples of conformal blocks for higher values of c from tau functions associated 
to branched covers of P 1 [24] . 

It is in principle straightforward to obtain from (1.8)— (1.11) similar expansions for 
Painleve V and Painleve III tau functions. In particular, this gives full short-distance 
(conformal perturbation theory) expansion of two-point functions in the Ising and free- 
fermion sine-Gordon field theory, and of the PV tau function describing universal part 
of the process of formation of Fisher-Hartwig singularities in the asymptotics of Toeplitz 
determinants [8]. This also seems to shed some light on recent results of [1, 17]. We hope 
to return to these issues elsewhere. 
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